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Graphene is promising as a host material for electron spin qubits because of its predicted potential 
for long coherence times. In armchair graphene nanoribbons (aGNRs) a small bandgap is opened, 
allowing for electrically gated quantum dots, and furthermore the valley degeneracy is lifted. The 
spin lifetime T\ is limited by spin relaxation, where the Zeeman energy is absorbed by lattice 
vibrations, mediated by spin-orbit and electron-phonon coupling. We have calculated Ti by treating 
all couplings analytically and find that T\ can be in the range of seconds for several reasons: (i) low 
phonon density of states away from Van Hove singularities; (ii) destructive interference between two 
relaxation mechanisms; (iii) Van Vleck cancelation at low magnetic fields; (iv) vanishing coupling 
to out-of-plane modes in lowest order due to the electronic structure of aGNRs. Owing to the 
vanishing nuclear spin of 12 C, Ti may be a good measure for overall coherence. These results and 
recent advances in the controlled production of graphene nanoribbons make this system interesting 
for spintronics applications. 

PACS numbers: 87.75.-d, 76.60.Es, 63.22.Rc 



I. INTRODUCTION 

Graphene has attracted intense scientific interest for its 
mechanical, electronic, and other properties. 1 4 Within 
the plane of its two-dimensional lattice it is extremely 
rigid while out-of-plane deformations are relatively soft 
due to the lack of a linear restoring force. 5 6 The ab- 
sence of a bandgap leads to a quasi relativistic behav- 
ior of the electrons that can be described by a Dirac- 
like Hamiltonian. 7 8 However, for typical semiconductor 
applications like transistors or spintronics devices, it is 
favorable to work with a bandgap. 9 12 Due to Klein's 
paradox, a bandgap is necessary to confine charge carri- 
ers electrostatically in graphene. 7 13 There are different 
situations that lead to a bandgap in graphene and some 
of them have already been studied in view of spintronics 
applications. 12 14 

Armchair graphene nanoribbons (aGNRs) can exhibit 
a bandgap and in addition allow for coupling of qubits 
in non- adjacent quantum dots (QDs). 15 17 Such a non- 
local coupling of qubits is ideal for fault-tolerant quan- 
tum computing and thus for scalability. 12 18 Over the 
past years, there has been substantial progress towards 
the goal of controlling the GNR edge termination within 
the production process and the controlled production of 
aGNRs might become feasible in the near future. 19 20 

Spintronics applications like the Loss-DiVincenzo 
quantum computer require spin coherence times much 
longer than typical operation times. 21 22 When the qubit 
is represented by the real electron spin, carbon materials 
are considered promising due to the small atomic spin- 
orbit coupling and weak interaction with nuclear spins in 
carbon. 12 23 24 While the curvature significantly enhances 
intrinsic spin-orbit coupling and hence spin relaxation 
in carbon nanotubes, this effect should not occur in flat 
graphene. 23 25 28 The natural abundance of 13 C, the only 
stable carbon isotope with a finite nuclear spin 1 = 1/2 
is only 1%. The concentration of nuclear spins can be 
further decreased by depleting this isotope. For magnetic 



fields above the lOmT-regime, flip-flop processes between 
nuclear spins and electronic spins become suppressed due 
to the different magnetic moments, /iB ^> Mnuc- We ex- 
pect that T 2 is dominated by Ti and that the spin re- 
laxation time is a good measure for overall coherence, 
Ti « T 2 /2. 

In this paper, we calculate the spin relaxation time T\ 
for electrons that are confined in an aGNR QD. The fi- 
nite width of the quasi one-dimensional aGNR leads to 
confinement in the transverse (x-) direction. As we will 
discuss in Sec. Ill, aGNRs of appropriate width have a 
bandgap. This allows to avoid Klein's paradox and con- 
fine electrons in the longitudinal (y-) direction by means 
of an electrostatic potential V(y). In a perpendicular 
magnetic field Be zi the two possible spin states of an 
electron inside the QD are split by the Zeeman energy 
gii^B = hu, where g = 2 is the electron ^-factor. Fig. 1 
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FIG. 1: (Color online) Sketch of the system and definition 
of the coordinate frame. The GNR has armchair termina- 
tions in the x-direction. The width of the sketched aGNR 
is characterized by m = 3 and /x = — 1 which leads to a 
bandgap that allows for electrostatic confinement in the y- 
direction. The potential V(y) defines the two barrier regions 
Bl, B2 (shaded) and the dot region D, that lies symmetri- 
cally between the barrier regions. The interatomic distance 
in graphene is a = 1.42 A. 
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shows a sketch of the system. 

Due to energy conservation, the Zeeman energy must 
be transferred to the lattice upon spin relaxation. For 
typical lab magnetic fields B < 20 T, the Zeeman energy 
corresponds to low-energy acoustic phonons at the center 
of the Brillouin zone. 29 We consider two cases separately: 
(i) free and (ii) fixed boundaries. The electron-phonon 
coupling Hepc comprises the deformation potential as 
well as the bond length change and couples in-plane vi- 
brational modes to the electronic state. By including 
the spin-orbit interaction Hsoi, the spin thus becomes 
connected to the vibrational state of the system. The 
coupling to the out-of-plane modes is considered, as well. 
Yet such a coupling either vanishes identically due to the 
electronic structure in aGNRs or appears only in higher 
order. 

This paper is organized as follows: In Sec. II, we 
present our model and in Sec. Ill, we recapitulate the 
bound states of aGNR QDs and explain the extended, 
quasi continuous states. Acoustic GNR phonons are 
shortly reviewed in Sec. IV. The effective spin-phonon 
coupling mechanisms that lead to T-f 1 via Fermi's golden 
rule are clarified in Sec. V. In Sec. VI, we comment on 
the actual evaluation of T\ . The results are presented in 
Sec. VII and discussed in Sec. VIII. 



states. 30 For this admixture, we consider both bound 
states confined inside the dot and extended, quasi con- 
tinuous states energetically above the confinement poten- 
tial. 

As will be discussed in more detail, Hsoi perturbs 
the electron-spin product states \k)\s) = \ks)(°\ where 
s =t,i- We denote the first order perturbed states by 
\ks). Finally, the electron-phonon coupling leads to fi- 
nite matrix elements (k H'HepcI^ t)- This allows us to 
use Fermi's golden rule in order to calculate the spin re- 
laxation rate 

T i _1 = ^Z)KH,n a>g + i|?/EPc|fet,n a , g )| 2 

a,q 

X States (fi^a.g) , (4) 

where Estates (^a,q) is the phonon density of states at 
the respective energy. The result is a function of three 
parameters: (i) length-to-width ratio (aspect ratio) L/W 
of the QD, (ii) potential depth AV of the QD, and (hi) 
perpendicular magnetic field B. We find that T\ can be 
as large as several seconds if Estates is small and the two 
mechanisms in Hepc interfere destructively. 

III. ELECTRONIC STATES 



II. MODEL 



The Hamiltonian of the system is 

T-i = ^elec + ^phon + %SOI + ^EPC , 



(i) 



where H e \ec and H p hon describe the unperturbed elec- 
tronic system and the unperturbed vibrational system, 
respectively. The spin-orbit interaction Hsoi leads to an 
admixture of opposite spin states such that the electron 
phonon coupling Hepc can induce a spin flip. Denot- 
ing the Fermi velocity by and the pseudospin by tr, 
the unperturbed electronic part of the system obeys the 
Hamiltonian 



K 



elec " 



-ihv 



Fl -cr x d x +cr y d y 



-V(y) (2) 



with eigenstates \k). The pure vibrational modes are 
described by 



H 



phon 



a,q 



h^a,q 



v a,q 



(3) 



where the summation runs over all phonon branches a 
and wavenumbers q. The angular frequency uo aA of a 
vibrational mode is implicitly determined by a and q and 
n aA is the occupation number operator. The eigenstates 
are the occupation number states \n a ^ q ). 

Since Hepc does not couple to the spin, the spin-orbit 
interaction Hsoi needs to be included in order to obtain 
a spin relaxing mechanism via admixture of electronic 



Due to the aGNR edges where the wavefunction van- 
ishes on both sublattices, electronic states in an aGNR 
have transverse wavenumbers 



n(n-fi/3)/W, 



(5) 



where n = 0,±1,±2, ... and W = (3m + fi)y/3 a is the 
ribbon width. 15 The width depends on m G N and 
fi G { — 1,0, +1}. The interatomic distance is a = 1.42 A. 
Due to Eq. (5) and E = ±hvFy/q* + & 2 , where vp is 
the Fermi velocity and k is the longitudinal electronic 
wavenumber, there is a bandgap ^ gap = ^vf\qo\- Since 
Eg&p = for // = 0, we assume /i = ±1 from now on. 
Note that /i is determined by the number of atoms across 
the GNR, Fig. 1 . Spinors with different transverse quan- 
tum number n are orthogonal such that we shall focus 
on the lowest transverse wavenumber with \q \ = tt/3W. 
The resulting gap £? gap = 2hvFir/3W allows to avoid 
Klein's paradox and confine charge carriers electrostati- 
cally in a finite square potential 12 



V{y) 



-{ 



: y G D (dot region), 

AV : y G Bl U B2 (barrier regions). 



(6) 



The barrier region Bl extends from the left end of the 
aGNR to y = and the barrier region B2 extends from 
y = L to the right end. The dot region D lies symmetri- 
cally between the barrier regions. The resulting potential 
landscape is shown in Fig. 2 (a) together with a bound 
state, which will be discussed in the next subsection. The 
length of the QD is denoted by L and assumed to be 
much smaller than the overall ribbon length, L <C £gnr- 



3 









/ Pn \ 

E 




















Egap — ^hvpqo 




Bl 




f D 


B2 



„ L 




AV [hv F q ] 

FIG. 2: (Color online) Electron states in an aGNR QD. (a) 
Sketch of a bound state and (b) QD bound state energy spec- 
trum given by roots of Eq. (12). (a) Due to armchair bound- 
aries, the minimum transverse wavenumber is qo = zLtt/3W 
for fi = =pl. As a consequence, the conduction band is sepa- 
rated from the valence band by a gap of E gap = 2^f|^o|- All 
energies shall be measured w.r.t. the middle of this bandgap 
inside the QD region. In the barrier regions, both bands are 
shifted by the barrier height AV. The resulting QD hosts at 
least one bound state. All bound states have the form given 
by Eq. ( 11 ) and decay exponentially for y — > zboo. The arrows 
underneath the Greek letters indicate the directed character 
of the according part of the wavefunction. The plotted prob- 
ability density \ip(y) | 2 belongs to the lowest bound state for 
L/W = 5 and AV = 1.8hvFqo> (b) Bound states exist for 
roots of Eq. (12) and can be plotted in a AV-E-p\ot. There 
is at least one bound state for all values of AV. As Al/ is 
increased, more bound states fit into the energy gap until 
the lowest state can leave the QD via valence states in the 
barrier regions. Notably, Eq. (12) has exactly one root for 
every value of E > TivyQo- We enumerate the bound states by 
j = 0, 1, 2, . . .. The circled position on the j = line marks 
the state plotted in panel (a). For the shown plot, the aspect 
ratio is L/W = 5. 



For concreteness, we assume an overall GNR length of 
Lgnr = 50W. The finite square potential needs to be 
considered in the electronic dispersion relation, which be- 
comes 

E = V(y)±hv F y/q* + k*. (7) 



Provided that the barrier height AV does not exceed a 
critical value 2hv^\qo \ + AVi, we can easily order bound 
states and extended states by their energies. The critical 
value and AVi will be explained in the next subsection 
- for now, we only assume that AV does not exceed it. 
Then, a state with energy E G [7wf|(7o|> ^f|#o| + AV] 
is bound since its longitudinal wavenumber k is real in 
the dot region and complex in the barrier regions, thus 
leading to an evanescent behavior. For E > hv^\qo\-\-AV, 
the longitudinal wavenumber is real in all regions. This 
leads to extended waves. Both bound and extended states 
contribute to the admixture mechanism and thus shall be 
discussed in more detail. 



A. Bound states 

To describe bound states in aGNRs, one can assume an 
infinite ribbon. 12 On one hand, Lqnr will always be finite 
in reality. On the other hand, bound states are mainly 
localized in the dot region < y < L and decay expo- 
nentially in the barrier regions, as shown in Fig. 2 (a). 
As mentioned above, we assume L <C Lgnr, such that 
the overall ribbon still appears approximately infinite for 
bound states. This allows us to follow the description 
with Lqnr — > oo for bound states. 12 

Accordingly, we denote the four component envelope 
wavefunction by 

* = &p,*™,-i>f\^i>F ) ) (8) 

and assume plane waves along the ribbon, ip^(x^y) = 
X^k( x ) e±iky ^ where 

= 4 +) (1, 0, 0)e^ x + &(+> (-* n>jfc , 1, 0, 0)e-^ x 
+c« (0, 0, -z n ^ l)e^ x + d« (0, 0, 1, z ntk )e-^ x (9) 
and 

X { ~1 = a£"> (z n , k , 1, 0, 0)e^ x + 6^(1, -z n , k , 0, O^"*"* 
+4r\o, 0, 1, -z^yi™* + S~\0, 0, z„, fc , l)e- i «« x . (10) 

With z n ^ = ±(#n + ik) I \[<£ti + an d longitudi- 
nal wavenumbers /cd = \J (E/hvp) 2 — q% (dot region), 
= k^/i = — (C^ — eAV) /Ttvp) 2 (barrier re- 

gions), bound states have the form 

[ <XnxkLe KBV : ye Bl, 

^ = \ PnX { nL elk " V +7nxi;2 D e-^ : y G D, (11) 
{ S nX £Le-^y- L ) : ye B2. 

The matching conditions at the interfaces Bl/D and 
D/B2 (that is, at y = 0,L) are discussed in Ref. [12 
and can be met for roots of the transcendental equation 

tan(fc D L) = k pKB _ (12) 
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FIG. 3: (Color online) Sketch of an extended state. The 
potential landscape, the aspect ratio L/W, and the barrier 
height AV are the same as in Fig. 2 (a) for bound states. 
The plotted probability density belongs to an extended state 
that is incident from the left as described by Eq. (14) and for 
which /ceb = 207t/Lgnr. The arrows underneath the Greek 
letters indicate the direction of propagation of the according 
part of the wavefunction. 



ergy is conserved, the wavenumber becomes 



&ED 




q 2 + fc| B + AV/hvp 



Q 2 



(13) 



in the dot region, uepenamg on tne 
state is incident from y < 0, leading to 

w { — ) p — ikED 



Depending on the sign of /ceb, the 
n y < 0, leading to 



<H /3"xS ED e lfcEDy + 7 n xS ED 

^xS"L eJfcEB(y " L) 

incident from y > L, which is described 



ye Bl, 
yeD, (14) 
ye B2, 



see Fig. 3 , or it is 

by 



i Jn ^n,kEU ^ i n ^n,k E p 

c ( + ) ih^(n,— T,\ , 



,(+) pik E B(y-L) I -v( _ ) p-ifcEB(y-i) 



2/eD, 

2/GB2. 
(15) 



For | g 1 = tt/3W and L/1V = 5, Fig. 2 (b) shows these 
roots as a function of the barrier height AV. There is 
a finite number of longitudinal excitations for any given 
AV. The different bound states can be enumerated by 
j = 0, 1, 2, ... and have distinct coloring in our figure. The 
j-th bound state has j nodes inside the dot region. For a 
given band, AV can be increased until the valence band 
reaches the energy of the lowest state which can now 
leave the QD via valence states in the barrier regions. 
Note that this occurs exactly when the argument on the 
l.h.s. of Eq. (12) equals a multiple of 7r. For the ground 
state, this means &d £ [0,7r/L] such that the maxi- 
mum ground state energy is £?o,max = ^fv?o + i 71 /^) 2 - 
States of higher energy belong to the j-th longitudinal 
excitation (j > 0) which begins at AVj = hvp(j7r/L). 
For AV < AVi, the ground state is the only bound 
state. This will be important for the evaluation of 7\, 
see Sees. VI and VII. 

The critical value for AV mentioned before is AV = 
2hv-p\qo \ + AVi. If the barrier height surpasses this value, 
the lowest state inside the QD can leave it via valence 
states in the barrier region. That is the state becomes 
extended thus affecting the ordering of bound and ex- 
tended states. Throughout this paper we assume that 
AV does not exceed this threshold such that the ground 
state lies on the j = band. 



B. Extended states 

We assume Lqnr = 501V for the overall length of 
the GNR such that possible wavenumbers are &eb = 
0, ±27t/Lgnr, • -j ±7r/a with lattice constant a. Since en- 



The matching conditions at y = 0, L can always be met. 
In contrast to bound states, extended states are propa- 
gating waves in the barrier regions. 



IV. ACOUSTIC GNR PHONONS 

The phonon energies we are interested in need to match 
the Zeeman splitting, Tiuj = g\i^B, where uo is the phonon 
frequency, g the electron g-factor, and /ib denotes Bohr's 
magneton. For typical lab magnetic fields B < 20 T, 
this implies low-energy acoustic phonons at the center of 
the Brillouin zone which can be modeled by continuum 
mechanics. 29 32 In this model, deformations are described 
by the displacement field u(r). While the components 
u xz and u yz of the strain tensor = (diUk + dkUi)/2 are 
known to vanish for thin plates in general, the monatomic 
thickness of graphene implies that u zz must vanish, as 
well. With Ui Z = 0, the elastic Lagrangian density of 
monolayer graphene is given by 29 33 34 



C = T-V 



P -2 
2 U 



;(a^) 2 



B 



M 2 

Uaa 



where A = d 2 + d 2 



and u 2 ik 



the sum convention with uu = 
• • has been 



u 2 xx + v 2 



^XX H" Uyy \ ^ ZZ vu ik — | ^ xy 

used, p is the mass density, and n is the bending rigidity. 
The bulk (B) and shear (/i) moduli can be expressed by 
Poisson's ratio a and Young's modulus £. The numeri- 
cal values of elastic and other constants we use are listed 
in Table I, Eq. (16) shows that in-plane vibrations uu 
decouple from out-of-plane vibrations u±. By assuming 
Ui(x, y) = f l {x) exp[i(qy — cut)] for a single mode and im- 
posing free or fixed conditions as discussed in Ref. [29J, 
we obtain the according phonon dispersions (Fig. 4) and 



5 



g>0 
'3 5 



Qo 



(a) 

_ -"-"05 






^^0^4,2 ^^^^ 












free boundaries 



20 
16 
12 

8 



to 




modes cj X2/ is related to the physical frequency by gj^ = 
uj^p/ ( E W. In the continuum model, all branches ex- 
tend to infinity but we are only interested in the range 
uj xy G [0, 5] since for a typical GNR width of W = 30 nm, 
&xy — 5 relates to a magnetic field of 20 T. The di- 
mensionless wavenumber q is obtained from the physical 
wavenumber via q = qW. For symmetry reasons ex- 
plained in Sec. VI, not all branches contribute to spin 
relaxation but only those which are explicitly labeled in 
Fig. 4. 

In the case of free boundaries, the branches ol^ 
and a 5 are relevant. While a 2 extends throughout the 
considered interval, a 5 only exists above Cu xy = 4.16, and 
«4 needs further discussion. Its minimum is 3.05 and 
occurs at a finite value % where the density of states 
has a Van Hove singularity. We consider two parts: for 
q < qo, we label the branch a 4jl and its label for q > q 
is o<4 5 2. The range of is Cu xy G [3.05,3.18] and a ^2 
extends from its minimum to infinity. 

For fixed boundaries, we only need to consider the 
branch aq. It extends from uj xy = 2.06 to infinity. We 
emphasize that for a typical GNR width of 30 nm, single 
phonon processes do not occur up to 8.25 T as there are 
no phonons below uj xy = 2.06 for fixed boundaries. 



FIG. 4: (Color online) Phonon dispersion for (a) free and (b) 
fixed boundaries. The dimensionless frequency Cu xy is con- 
nected to the physical frequency by uj xy = uj^J p/*EW, where 
the radicand contains elastic constants listed in Table I We 
restrict our interest to the frequency range uj xy G [0,5] since 
for W — 30 nm, the upper bound already relates to a mag- 
netic field of 20 T. The scale on the r.h.s. shows the magnetic 
field for W = 30 nm. Due to parity w.r.t. x, only the labeled 
branches (a* for free and cti for fixed boundaries) assist in 
spin relaxation, (a) The phonon spectrum is gapless for free 
boundaries. The branch ol± has a minimum and hence a di- 
verging density of states for finite q. Its constituent parts a^i 
and «4,2 shall be treated separately, (b) Fixed edges lead to 
gapped phonon spectrum. For W = 30 nm, this gap corre- 
sponds to 8.25 T. In our range of interest, the branch labeled 
ol\ provides the only channel for spin relaxation. 



the explicit displacement fields. The latter can be quan- 
tized and then take the form 



t x P + fv P \ P i( iy 
u ± =J2r a , q fa, q e z e^, 



/ j ' Oi,qyj a, 1 
a,q 



(17) 
(18) 



a,q 



where q is the phonon wavenumber, a labels the phonon 
branch, and 



h/(2pLWu; a , q )(b a , q + b^_ q ) 



(19) 



is the normal coordinate. The operator b a ^ q (fot ) annihi- 
lates (creates) a phonon on branch a with wavenumber q. 

As discussed in the following section, we can ne- 
glect coupling to out-of-plane modes and focus on in- 
plane modes. The dimensionless frequency of in-plane 



V. COUPLING MECHANISMS 

Phonons do not couple to the electron spin directly. 
The relevant mechanism usually involves the spin-orbit 
interaction. 23 24 30 In graphene, the spin-orbit interaction 
is given by 



H^oi — Hi + Hr 

= ^IT Z (J Z S Z + \R(T z (J x S y - (7 y S x ) 



(20) 



and we will consider it in order to obtain an indirect spin- 
phonon coupling. 26 28 The intrinsic (Dresselhaus) term 
Hi has coupling strength A/ and the Rashba (or extrin- 
sic) term Hr couples with strength Xr. The valley is 
denoted by t z , pseudospin by cr, and real spin by s. In 
the following, we show how the real electron spin can 
be connected to the vibrational state of the system by 
taking the spin-orbit interaction into account. 



A. Coupling to in-plane modes 

In first order perturbation theory, Hr corrects the 
electron-spin product states \k)\ J [) = \k^)^ to 



ifct) = ifct) 



(0) 



fej ki) E k -E> k+ g» B B> (21) 



and I k \) accordingly. We emphasize that the summation 
index k' runs over both bound and extended states. The 
potential depth and the aspect ratio determine how many 
bound states exist, Fig. 2 (b). For extended states, we 
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consider all wavenumbers inside the first Brillouin zone, 
&eb = 0, ±27t/L G nr, ±7r/a. 

The second term in Eq. (21) admixes states with op- 
posite spin such that the electron-phonon coupling Hepc 
can induce a spin flip 30 



(fe||«EPc|fet> 

\Umc)kk>{U R )$ k 



(H EP c)k>k(H R )H 
E k - E k > + gusB ' E k - E k , - g/j, B B 



+ 



(22) 



where we denote the numerator in Eq. (21) as (Hr)^}, 



k'k 

and the spin-conserving transitions of Hepc accordingly. 
We find that for a given k' the two terms in Eq. (22) ex- 
actly cancel each other at B = 0. This effect is known as 
Van Vleck cancelation and is expected for time-reversal- 
symmetric systems. 35 Moreover, (Hr)]^., vanishes if both 
k and k' represent bound states and the longitudinal ex- 
citation indices j^, jV (see Fig. 2) are both even or both 
odd. In the electron phonon coupling Hamiltonian Hepc, 
we consider the deformation potential Hdp as well as 
bond length change Hblc' 



Hepc = ^dp + %blc , 
%dp = 0i V • U\\ , 

/ A 

%BLC = 92 



(23) 




where gi^ are coupling constants, (A x ,A y ) = (u xx — 
u yy , —2u xy ), and the basis of Eq. (8) has been used. 8 33 34 



B. Vanishing out-of-plane deflection coupling 

Low-energy acoustic phonons at the center of the Bril- 
louin zone have a wavelength much larger than the lattice 
constant and produce a local tilt of the GNR. In the lo- 
cal ribbon frame E' where n = e' z is the vector normal 
to the ribbon plane the local spin matrix is described 
by s' z = s z — d x u z s x — d y u z s y . As a consequence, the 
intrinsic spin-orbit interaction 



Hi = \iT z a z (s z - d x u z s x - d y u z s y ) 



(24) 



becomes dependent on out-of-plane phonons such that 
these could flip the spin. This is known as deflection 
coupling. 24 31 However, there is a proportionality to r z 
in Eq. (24). Since the electronic states we use have the 
property that the wavefunction has equal weight on each 
valley, 36 the contributions from K and K' add up to zero 
and the deflection coupling between spin and out-of-plane 
modes vanishes. 

If the spin-orbit admixed states of Eq. (21) are used, 
there is a finite overlap only between both admixed parts 
such that the resulting mechanism is proportional to 
A/A 2 ^ and hence negligible. 



Compared to in-plane phonons, both deformation po- 
tential and bond length change appear only in higher 
order such that we neglect these mechanisms for out-of- 
plane phonons. 



VI. EVALUATION OF Ti 

Using Eq. (4), we calculate the spin relaxation rate for 
the electron in the lowest bound state (ground state) of 
the QD. For concreteness, we assume /i = — 1 in Eq. (5). 
According to Eqs. (11), (14), and (15), both bound and 
extended states have a finite probability density in the 
barrier regions. However, bound states are localized in 
the dot region and decay exponentially in the barrier re- 
gions. In particular the ground state, plotted in Fig. 2 
(a), has a very low probability density in the barrier re- 
gions. Its overlap with another bound state in the barrier 
regions is negligible. Only high-energy extended states 
have a significant contribution outside the dot region. 
Yet the overlap of an extended state with the ground 
state outside the dot region will still be small and since 
they are energetically far apart, the contribution from 
the barrier regions can be neglected. As a consequence, 
we can restrict the integrals in Eq. (22) to the dot region. 

As discussed in Sec. IV, we consider phonons with free 
boundaries as well as phonons with fixed boundaries. Not 
all branches contribute to spin relaxation: Because of 
mirror symmetry w.r.t. x = 0, Hblc and Hbp are even 
or odd in x, depending on what phonon branch they be- 
long to. Due to their similar form, 37 the mechanisms are 
either both even or both odd for a given branch. The 
^-dependencies of the electronic states in the matrix ele- 
ment (T^EPcOfc'fc cancel out, e l< ^ qo ~ q ^ x = 1, such that the 
x-integral vanishes if Eq. (23) is odd in x. The branches 
«2, «4, and a$ in Fig. 4 (a) and a\ in Fig. 4 (b) have 
couplings Hepc that are even in x and hence can relax 
the spin. 

For a given relaxation channel (a,q), both mecha- 
nisms Hdp and Hblc are combined in Eq. (4) coherently. 
Moreover, the couplings via bound states and extended 
states in Eq. (21) are added up in a coherent way. We 
are interested in the relaxation of the spin in the ground 
state, which corresponds to j = in Fig. 2 (b) and hence 
restrict the barrier height to AV G [0, 2hvFqo + AVi]. 
If AV exceeds this upper bound, valence states become 
available in the barrier regions and the lowest state in- 
side the QD can leave the dot region. For AV < AVi on 
the other hand, the ground state is the only bound state 
such that the perturbation in Eq. (21) comes about only 
due to extended states, which fully determine the spin 
relaxation in this case. 

For spin relaxation, Eq. (4) is proportional to n a ^ q + 1 
and we assume n aA = 0, i. e. sufficiently low tempera- 
ture, k^T <^huj = qhbB' By &b we denote Boltzmann's 
constant. Assuming a magnetic field of B = 1 T, this 
means T <C 1.3 K. For T > 15 K, spontaneous emission 
can be neglected since n aA ^> 1 and one obtains the spin 



relaxation by multiplying our results with the expecta- 
tion value of the Bose distribution 



B(w=: 



(VK=30nm) 



[T] 



(n a , g (B,T)) 



(25) 



The spin relaxation time T\ is a good measure for over- 
all coherence when pure dephasing, which comes from 
coupling to nuclear spins, is negligible. Due to the low 
density of nuclear spins in natural carbon and the very 
different magnetic moments [1b ^> /inuc? we expect that 
flip-flop processes with nuclear spins can be neglected for 
magnetic fields above 10 mT. For a typical GNR width 
of W = 30 nm, 10 mT correspond to uo xy = 0.0025. As 
a consequence, we restrict our calculations to the inter- 
val Cu xy G [0.0025,5]. The upper bound corresponds to a 
magnetic field of 20 T. All plots that show rates are cut 
off at these boundaries. 



VII. RESULTS 

To calculate Ti, we need to use the specific values 
of the elastic constants that define the phonon spec- 
trum. Young's modulus for the two-dimensional lattice 
of graphene is obtained by multiplying the bulk value 
with the thickness associated with graphene, <E = £3D^> 
where h = 3.4 A. For further discussion of the elastic con- 
stants, we refer to Ref. [29 . Table I gives an overview 
of the constants we use in our calculation. The Rashba 
type spin-orbit coupling is linear in the electric field and 
thus can be adjusted by an external electric field or by 
using a suitable substrate. 27 



a = 0.16 
-E 3.4TPaA 
B = 12.6 eV/A 2 
/i = 9.1eV/A 2 



38-40 gi — 30 eV 24 33 34 

38 40 #2 = 1.5 eV 24 34 

6 40 A R = 40 x 10" 6 eV 24 27 28 

6 40 vf — 8.8 x 10 5 m/s 8 24 28 



7.61 x 10~ 7 kg/rn 4i 



TABLE I: Numerical values of the parameters we use in our 
calculation. 

The spin relaxation time T\ depends on three parame- 
ters: (i) the aspect ratio L/W of the QD, (ii) the poten- 
tial depth AV of the QD, and (hi) the applied perpen- 
dicular magnetic field B oc u xy . Moreover, the phonon 
spectrum and hence the spin relaxation depends on the 
mechanical boundary conditions. We discuss free bound- 
ary conditions separately from fixed boundaries. 




FIG. 5: (Color online) Partial rates for various relaxation 
channels. For L/W = 5 and AV = l.ShvFqo, all contributions 
to the four relaxation channels as (a), 0:4,1 (b), 04,2 (c), and 
«2 (d) are shown. The contributions stem from 7ir>p with 
admixture of bound (labeled "bDP" ) states or extended states 
("eDP") and from 'Hblc with the same admixtures ("bBLC" 
and "eBLC", respectively). These contributions are added 
up coherently to the "combined" relaxation of the respective 
channel. 



A. Free boundary conditions 

For symmetry reasons explained above, only the 
phonon branches with labels 0:2, <^4 (consisting of parts 
c^4 5 i and 0^4,2)7 and a$ in Fig. (4) (a) need to be consid- 
ered. The respective rates of these relaxation channels 



are shown in Fig. 5 for an aspect ratio of L/W = 5 and 
a barrier height of AV = LShv^qo. The allocation of 
panels to branches is as follows: panel (a) belongs to 
branch a 5 , (b) to o^i, (c) to <%4 5 2, and (d) to oli> Each 
panel shows four separate contributions to T x _1 that come 
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about from the two mechanisms in Eq. (23) and the ad- 
mixture of bound states or of extended states in Eq. (21) 
for each mechanism. The coherent sum of all four con- 
tributions is displayed by the grey line. The deformation 
potential usually dominates over the bond length change 
since its coupling constant is 20 times larger, Table I, 
For AV = 1.87m;f(7o> extended states are energetically 
far away from the ground state such that the contribu- 
tion from the deformation potential with admixture of 
bound states dominates in Fig. 5. Oscillations in indi- 
vidual rates may be due to the phonon phase e iqy that 
is integrated with the matrix elements (Hepc) k'k and 
rotates according to the phonon dispersion when lj is 
changed. Panels (c) and (d) show that the matrix el- 
ements (Hr>p)k'k an d (T-Lblc) k'k ma y interfere destruc- 
tively, thus decreasing T x _1 by several orders of magni- 
tude, yet typically not to zero. 

In all these plots, the bottom scale shows uj xy and the 
top scale shows the magnetic field B that corresponds to 
uu xy , assuming a width of W = 30 nm. 

Fig. 6 (a) shows the full spin relaxation rate for the sit- 
uation of Fig. 5, that is, the combined rates of all relax- 
ation channels o^, <^4,i, <^4,2> and a$ (grey lines in Fig. 5) 
are summed up to the full relaxation rate T-f 1 (grey line 
in Fig. 6 (a)). The rate with the label "bound" ("ex- 
tended") is obtained in a similar fashion, but only con- 
tributions with admixture of bound (extended) states are 
considered, here. For AV = l.Shvpqo, the admixture of 
bound states dominates the admixture of extended states 
by several orders of magnitude. Yet by lowering AV, 
the influence of extended states can be close to (Fig. 6 
(b)) or even surpass the influence of the bound states. 
Fig. 7 shows two cases for an aspect ratio of L/W — 2. 
In panel (a), the barrier height is AV = 2.0^vf<7o and 
extended states are basically irrelevant compared to the 
relaxation via bound states. However, panel (b) shows 
that for AV = 0.9hvpqo, the major contribution comes 
from the extended states. 

Fig. 8 (a) shows T-f 1 as a function of parameters AV 
and O xy oc and for a fixed aspect ratio of L/W = 5. 
In contrast to Q xy , the barrier height hardly changes the 
qualitative picture. The orange cut at AV = l.Shv-pqo 
is repeated in panel (b) in a doubly logarithmic plot 
that highlights the B 5 dependence in the range uj xy £ 
[0.0025,0.5]. In this range, only the branch ot2 is avail- 
able and has a linear dispersion B oc uj oc q. The matrix 
elements (Hepc) k'k have one power in B due to (i) the 
gradients cx q in Eq. (23), (ii) dipole approximation cx 
and (iii) Van Vleck cancelation cx B, each. Because of the 
prefactor cx cj" 5 in Eq. (19), we find (%EPc)fc'fc oc B 2 5 . 
As OL2 is linear and hence p s tates °t B° for this low-field 
regime, this explains T-f 1 cx B 5 . Destructive interference 
of matrix elements {T-Lp>p)k'k an d (T-Lblc) k'k can lead to 
a very small but nonzero relaxation rate. 

Fig. 9 (a) shows a plot similar to Fig. 8 (a), yet for 
L/W = 2. The qualitative picture is much different from 
the aspect ratio L/W = 5. Figs. 6 - 9 (a) all show discon- 
tinuities at uj xy = 3.05 that stem from the branch for 



£>(W=30 nm) [T] 

4J) 8J) 12J) 16J) 20.0 




FIG. 6: (Color online) The relaxation rates for different dot 
depths AV. By summing up the combined relaxation rates 
(see Fig. 5) of all channels available for a certain 0J xy , the 
full relaxation rate (grey line) is obtained. The lines labeled 
"bound" and "extended" are obtained in a similar way by 
considering only bound or extended states, respectively. At 
& X y = 3.05, Tf 1 is discontinuous due to the advent of the 
relaxation channel that has a diverging density of states 
at this point, Fig. 4 (a). Panel (a) accords to parameters 
L/W — 5 and AV = 1.87m;f<?o as in Fig. 5 Clearly, the 
energetically far off extended states play a negligible role for 
such a deep dot. In panel (b), the barrier height is reduced 
to 0.27m;f<7o such that extended states are about as important 
as bound states. 



which the density of states has a Van Hove singularity at 
<2o while the coupling Hepc remains finite, Fig. 4 (a). 



B. Fixed boundary conditions 

Most importantly, fixed boundaries result in a gapped 
phonon spectrum. This means that spin relaxation in- 
volving only one phonon cannot occur for magnetic fields 
that correspond to uu xy < 2.06. Note, that for a typical 
width W = 30 nm, uj xy = 2.06 corresponds to a mag- 
netic field of 8.25 T. However, phonon scattering may 
still take place below this threshold. In contrast to our 
claim in Ref. [29 , only the branch 6l\ contributes to the 
spin relaxation rate. Gradients, dipole approximation, 
and Van Vleck cancelation play the same role as for free 
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£>(W=3Q nm) [t; 

4J> 8,0 12J) 16J) 20.0 




FIG. 7: (Color online) This plot shows the same quantities 
as Fig. 6, yet for the aspect ratio L/W — 2. Again, the 
influence of extended states depends on the barrier height: 
AV = 2.0hv F qo in (a) and AV = 0.9hv F qo in (b). The ex- 
tended states dominate in the latter case. 

boundaries, yet due to the gap the frequency uj is not pro- 
portional to some power of q such that there is no power 
law that connects T x _1 and B as for free boundaries. 

Fig. 9 (b) shows an analog to panel (a), yet for fixed 
boundaries. For aspect ratios larger than in Fig. 9 (b), 
oscillations occur which can again be explained with the 
phonon phase e iqy that rotates according to the phonon 
dispersion when uj changes. These oscillations arise only 
if the dot length is large enough. 



VIII. DISCUSSION 

The spin relaxation times we find in our work range 
from 10 _T seconds to beyond the second range. For 
cases where T\ is very long, it can be expected that other 
mechanisms not considered here will dominate. Our re- 
sults depend on the aspect ratio L/W, the barrier height 
AV, and the Zeeman splitting g^B ex uj xy but also on 
the mechanic boundary conditions that lead to different 
phonon dispersions. By choosing/adjusting these degrees 
of freedom properly, T\ can be in the range of seconds. 
We attribute such long relaxation times to several effects: 

(i) GNRs are quasi one-dimensional systems similar 
to carbon nanotubes. Both the phonon and the elec- 




0.1 0.2 0.5 1.0 2.0 5.0 



FIG. 8: (Color online) Spin relaxation rate T^ 1 for an aGNR 
with aspect ratio L/W — 5 and free edges, (a) The rate is 
shown as a function of barrier height AV and phonon fre- 
quency uj xy . The orange cut corresponds to the grey line in 
Fig. 6 (a) and is repeated in panel (b) with a doubly logarith- 
mic scale that highlights the B 5 dependence in the interval 
Co X y e [0.0025,0.5]. 



tron density of states are thus limited compared to bulk 
graphene. 23 

(ii) Destructive interference between the deformation 
potential and the bond length change as well as oscilla- 
tions due to the phonon phase e iqy that rotates according 
to the phonon dispersion when uj changes both reduce 
the relaxation rate T-f 1 by several orders of magnitude 
for specific magnetic fields. 

(iii) In contrast to other graphene QD systems, 
the electronic states in aGNRs are invariant under 
time reversal symmetry which leads to Van Vleck 
cancelation. 12 24 35 As a result, Eq. (22) vanishes for 
B = 0. 

(iv) Deflection coupling to out-of-plane modes vanishes 
as the evenly distributed weights on K and K' spinor 
components cancel out. As a result, only the very rigid 
in-plane modes need to be considered. This rigidity leads 
to a generally small density of phonon states p s tates- 12 31 
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FIG. 9: (Color online) The relaxation rate T^ 1 for (a) free 
and (b) fixed mechanic boundaries, (a) This case is similar to 
Fig. 8 (a) yet with aspect ratio L/W = 2. (b) Fixed boundary 
conditions and L/W = 2. Due to the gapped phonon spec- 
trum, the rate T-f 1 vanishes below uj xy — 2.06 for our model 
and with fixed boundaries. Moreover, T-f 1 is not discontin- 
uous in uj xy as the branch ol\ never becomes flat for finite q, 
see Fig. 4 (b). 



(v) Phonons do not couple to spin directly so that spin- 



orbit coupling needs to be included. However, spin-orbit 
coupling in graphene is weak compared to other systems 
(e.g. carbon nanotubes). 12 23 28 

(vi) The admixture of electronic states in Eq. (21) in- 
cludes bound and extended states. However, only every 
second bound state contributes; for parity in ^-direction, 
jk and jk' may not be even or odd at the same time, 
Fig. 2 (b). States that are energetically far apart from 
the ground state play a small role in the sum which is usu- 
ally the case for extended states, depending on AV. As 
a consequence, the admixture of these electronic states is 
suppressed. 12 30 

(vii) Due to parity in x-direction, not all phonon 
branches contribute to spin relaxation but only those 
with explicit labels in Fig. 4, for which Eq. (23) is even 
in x. This limits the number of relaxation channels. 29 
It is an open question how strong the avoided relaxation 
channels contribute to T x _1 if this symmetry is broken. 

(viii) We assume phonon vacuum in Eq. (4). A fi- 
nite temperature can be included by multiplying the rate 
T-f 1 with the expectation value of the Bose distribution 
(n a , q (B : T)} as explained in Sec. VI. 

The carbon isotope 12 C has no nuclear spin and the 
natural abundance of 13 C, which has spin 1/2, is only 
1%. Thus, pure dephasing which comes from coupling to 
nuclear spins is likely to play a minor role in graphene de- 
vices and Ti ~ T2/2 becomes a good measure for overall 
coherence. Our results show that electronic spin qubits 
in aGNRs are promising for spintronics applications like 
the Loss-DiVincenzo quantum computer. With view to 
recent advances in controlling the edge termination of 
GNRs it will be interesting to see whether aGNR spin- 
tronics can be realized in experiment. 
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